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^ \ We shortly review the progress in the domain of deterministic chaos for 

' quantum dynamical systems. With the appropriately extended definition of 

quantum Lyapunov exponent we analyze various quantum dynamical maps. 
. It is argued that, within Quantum Mechanics, irregular evolution for prop- 

erly chosen observables can coexist with regular and predictable evolution of 
, states. 
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The theory of deterministic chaos has a venerable history, dating back to Poincare cele- 
Qh| brated work (see ||50|| ). Here, deterministic chaos refers to deterministic development with 
chaotic outcome. Another way to say this is that from moment to moment the system is 
evolving in a deterministic way, i.e., the current state of a system depends on the previous 
state in a rigidly determined way. However, measurements made on the system do not al- 
low the prediction of the state of the system even moderately far into the future. It was 
.£h I observed that whenever dynamical chaos was found, it was accompanied by nonlinearity 
(cf. |P0|| ). Another observation made long time ago was that exponential sensibility in non- 
linear classical systems can lead to exceedingly complicated dynamical behavior (cf. |[37|| ). 
The characteristic features of such a behavior, e.g. deterministic unpredictability, positive 
Lyapunov exponents, are typical attributes of chaotic (deterministic) dynamical systems 
(cf. |]30|| , [J25|J ). We emphasize that the theory of such classical dynamical systems is a well 
developed subject and contains many (concrete) examples of chaotic systems. For instance, 
we can describe the dynamics of a classical system as the evolution in the phase space gov- 
erned by the canonical equations of motion. Solutions of these equations define the family 
of trajectories on the phase space. The complex behavior of certain class of trajectories 
can be reflected by the positiveness of Lyapunov exponents and for some choices of the 
Hamiltonians one can find concrete examples of such chaotic systems (see [|30| ). 

However, contemporary science has been founded on quantum mechanics. Therefore, it is 
quite natural to look for the quantum-mechanical description of (classically) chaotic systems 
(see e.g. 0, JTTJ] , 0, |0|, @, @, @, [0], Q). To achieve this objective it is 



important to realize that the chaotic behavior arises from nonlinearities in the equations of 
motion. Let us remind that the Schrodinger equation of motion is a linear one. Furthermore, 
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the Koopman's construction, in the theory of abstract dynamical systems, leads to linear 
evolution operators (see e.g. ||). On that basis, the long debate on the existence of quantum 
chaos leads to a conjecture that quantum mechanics suppresses chaotic behavior in such a 



way that there is only a room for "fingerprints" of chaos [12[ . 

Our objective in this paper is to show that quantum mechanics admits, in some special 
cases, certain features of chaotic behavior. To understand the origin of this phenomenon let 
us remind that the dynamics of a classical system can be described either as the evolution 
in the phase-space by the conventional Hamilton's equations or by the Liouville equations. 
Now if we consider a system with some ergodic properties e.g. mixing, then we can find 
that the chaotic and unpredictable behavior of some observables coexists with a regular and 
predictable behavior of densities (so states). The time evolution of a quantum system is 
described either by Schrodinger equations (quantum counterpart of the Liouville picture) 
or by the Heisenberg equations (quantum counterpart of the Hamilton's picture). We shall 
indicate how these scheme may be used to explain the just mentioned coexistence now in 
the quantum mechanical description. In other words, we shall show that a kind of un- 
predictability can also exist for quantum maps. To study this phenomenon we shall use 
quantum Lyapunov exponents. To make this point comprehensible let us point out that we 
shall study the stability of dynamical maps without any analysis of the integrability of the 
corresponding equations of motion. 

In Sections II and III we review the structure of classical mechanics and basic properties 
of (classical) chaotic systems respectively. This material is mostly well known; however, we 
emphasize the role of the picture in the description of the time development of systems. 
Section IV deals with the question of (non)linear lifting of dynamical maps. The objective 
here is to show the possibility of nonlinear lifting of such maps as well as its relation to cer- 
tain algebraic structure. Section V provides a detailed exposition of quantum counterparts 
of the previously reviewed characteristic features of classical dynamical systems. Section 
VI contains the definition of quantum characteristic exponents and their basic properties 
while Section VII indicates how these exponents may be used to study stability of quantum 
dynamical systems. We close, Section VIII, with a brief discussion on the obtained new re- 
sults. For the reader's convenience we include the Appendix with some basic " C*-algebraic 
vocabulary" . 



II. THE STRUCTURE OF CLASSICAL MECHANICS. 

We start with a review of the essential facts of classical mechanics. Let T = 
{q±, ...,q n ,pi, ...,p n } be the phase space of some classical system. The Newtonian systems 
with (local) forces derivable from a potential form the important and large class of dynamical 
systems. The dynamics of this class is described by the Hamilton's equations: 

. dH dH 

> Pi = --K— i = l,...,n (1) 
dpi dqi 

where qi,Pi denote generalized coordinate and momentum respectively, and H = 
H(qi, ...,q n ,pi, ...,p n ) is the Hamiltonian of the system. Let Q C T be a subset. We denote 
by A the set of all complex-valued smooth functions on Q. A is the set of classical observ- 
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ables of the considered system. The Hamilton's equations (||) lead to one parameter group 
of maps of the phase space into itself 



G t ■. r - r (2) 

where Gtfi, n € T, is the solution of the Hamilton's equations with the initial condition 
G t fi\ t=0 = Hq. The transformations G t can be lifted on the level of observables in the 
following way 



U t : A -> A (3) 



where 



Utffr) = f t (fi) = f(G tf x) (4) 

for any \i E V. Obviously, Ut is the one parameter family of linear maps on A. The prescrip- 
tion given by (01) is sometimes called the Koopman's construction. The next observation 
is 

gj&O _ y dft OH df t dH _ 

Here {■, ■} denotes the Poisson bracket which can be also considered as a derivation on the 
abelian algebra A of classical observables (cf. Appendix). In other words, if we equip A 
with the structure of abelian C*-algebra then the Poisson bracket determines one parameter 
group of automorphism of A and this group of automorphism represents the Hamilton's 
evolution of observables. Another basic concept in the description of classical systems is 
that of a state. Namely, a state of a system is described by a (probability) distribution g on 
T. In particular, a pure state is given by 

g(qi,...,q n ,Pi, -,p n ) = %i - Qi) ■ ••• • S(q n - gjj) • 5(p 1 -p?) • ... • S(p n -p°) (6) 

where 8(x) is the Dirac's 5-function, ((ft, q n ,Pi, ■••,p n ) is an arbitrary point in T while 
the considered fixed point in T is denoted by (g°, g°,p5, ...,p°). We remind the reader 
that evolution in classical mechanics can be given in one of two pictures. The one just 
described is so called the Hamilton's picture, i.e. time evolution of an observable is given 
by (0) while a state is time independent: 4f = 0. The second one, so called the Liouville 
picture determines the evolution of a state, i.e. the time development of a distribution g. 
The relevant equations of motion are 

Obviously, on the set {g — an arbitrary distribution} = S of distributions (states) one can 
define maps V t (in the same way as in (f|)) 

V t Q = g t - (8) 

Clearly, V t (agi + (1 — a)g 2 ) = aV t (gi) + (1 — «)^(^ 2 ) for any a e (0, 1). We remind that 
these pictures are equivalent in the sense 
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< ft > g = J r f t {fi)Q(jj)dfi = y r f{n)ot{n)dn =< f > et . (9) 

Thus, the following picture is emerging. The time evolution in classical mechanics can be 
described by either linear (U t ) or convex (Vt) maps. Therefore, there does not exist a room 
for nonlinearity which is a necessary condition for a chaotic behavior. But we remind (cf. 
Introduction) that there are many examples of concrete models with chaotic behavior. To 
explain this apparent contradiction let us make an important observation. Although the 
Hamilton's and Liouville pictures are equivalent in the sense given by (||) the Hamilton's 
picture offers larger possibilities. Namely, one can pick up a coordinate (or momentum) as an 
observable. Then the equations of motion given by (0) are reducing to that of the form (ffj). 
It is important to realize that equations ([!]) are on the phase space T (or on its subset) and 
for some choice of the Hamiltonian H(q,p) one can obtain a system of nonlinear differential 
equations. To understand this observation within the above presented structure of classical 
mechanics it is enough to note that to get the mentioned nonlinear equations we did not 
use the linear lifting, i.e. the Koopman's construction. Therefore, for a particular choice 
of (classical) observables and the Hamiltonian, there exists a chance to get a nonlinearity, 
so also to get trajectories with unpredictable behavior. There is another possibility to get a 
nonlinearity for classical systems. Namely, one can ask for a non-linear version of Koopman's 
construction. However, this is a more subtle point and we will discuss this question later on. 



III. CHAOTIC (CLASSICAL) SYSTEMS. 

Let us consider a system of differential equations given by (H) where for simplicity we 
put n = 1. Suppose R9tw x l (t) = (q l (t),p l (t)) E V, i — 1,2 is a solution of ([!]) with 
initial conditions ^(0) = (<f (0), p*(0)), i = 1,2. Put Ax(t) = dist(x x (t) , x 2 (t)) . Then one 
can ask the following question: what is behavior of the function IR E t i— > Ax(t)7 We are 
interested in the class of systems such that the growth of Ax(t) is of the exponential type. 
We remind that the exponential growth is measured by Lyapunov exponents which can be 
defined as the following limits: 

i) for a discrete dynamical system: 

X cl (x,y)= lim -\og\D x r n (y)\ (10) 



n.^oo 



n 



where D x r n (y) denote the directional derivatives of r composed with itself n times at a 
point x in a direction y (cf. p5|). 



ii) for a continuous dynamical system: 

X c \x,y)= lim J log \D x r t (y) | (11) 

where D x r t (y) denotes, as above, the corresponding directional derivative. 

Remark 1 If we examine a time evolution r t of a point x depending on a parameter v, then 

limt^oa^log\-^-T t (x(v))\ (12) 
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also can be considered as a measure of sensitivity of evolution with respect to initial condi- 
tions (determined by v) provided that the limit in exists. This limit can be also called 
Lyapunov exponent. 

Both limits measure the average rate of exponential growth of separation of orbits which at 
time zero differ by a small vector. This property is used as a measure of the sensitivity of 
dynamical system with respect to initial conditions. 

Definition 1 A dynamical system described by differential equations of the type ([][) is called 
regular (irregular) if X cl < (X cl > respectively). 

Lyapunov exponents are closely related with Kolmogorov-Sinai entropy £r-s which is an- 
other key quantity for analysis of classical dynamical systems. £r-s measures mixing prop- 
erties of a system and its relation to Lyapunov exponents for smooth enough systems is 
givenby (cf. [0, H, H) 



s K _s= E x * ( 13 ) 

k;X k >0 

Thus, for a large class of dynamical systems, a positivity of Kolmogorov-Sinai entropy indi- 
cates the existence of positive Lyapunov exponents. We recall that for more general class of 
dynamical systems (that is satisfying weaker smooth conditions) one should replace equality 
in (p~3|) by an inequality (see |36| ). 



Example 1 Let us consider a dynamical system determined by the following differential 
equation 

d 2 

—x = k 2 x, K>0. (14) 

til/ 

The solution of ftTft ) can be written as x(t) = Xicosh(nt) + X2sinh{nt) where xi, x 2 are 
constants describing the initial conditions. Put, for simplicity, x 2 = and assume that we 
have two initial conditions x[ and x[ J such that \x{ — x[ J \ = e. The corresponding solutions 
will be denoted by x\t) and x n (t) respectively. We note 

Ax(0) = 1^(0) - x H (0)\ = \x[ - xi 1 1 = e (15) 

while 

Ax(t) = e cosh(nt) - e Kt (16) 

Thus, the considered system exhibits the exponential growth of separation of trajectories, i.e. 
a sensitivity with respect to initial conditions. However, it is important to realize that it is 
difficult to consider this system as a chaotic one. 

This example clearly indicates that a definition of chaotic system presents a delicate problem. 
In particular, in the above example, we have gotten a sensitivity with respect to initial 
conditions since we allow trajectories to spread to infinity. This motivates 

Definition 2 A (classical) dynamical system exhibits a chaotic behavior if 
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• it is defined on a compact subset of the phase space (or the algebra of observables is 
defined on a compact subset of the phase space), 

• it has an irregular motion with some ergodic properties. 

Let us comment on this definition. Firstly, the assumption of compactness excludes a possi- 
bility of a system with trajectories going to infinity. In other words, an irregular evolution 
should have a kind of repeatibility. Let us emphasize that all computer simulations of chaotic 
systems have this feature. Next, we wish to briefly comment on the assumption of ergodicity. 
In real dynamical systems there are coexisting regions of chaotic behavior with regions of 
regular motion. So we need a kind of ergodicity to get a system with overwhelming region 
of irregular motion. 

We want to close this Section with the two remarks on the just given definition of 
chaotic system. Firstly, we want to point out that the discussed definition is not the only 
one possible. On the contrary, there are many different approaches to this topic with various 
definitions of chaos. Even for semigroups of linear operators on a Banach space, there is an 
attempt to introduce a concept of chaos (see the definition given by Deveney ||22|| , see also 



However, we strongly believe that for real mechanical systems the presented definition 
is the proper choice as it reflects our intuition on sensitivity of dynamics. Secondly, we 
want to stress again that the Liouville picture seems to be not the best choice for doing a 
study of chaotic systems. Clearly, one can examine the basic features of the Peron-Frobenius 
operator. Nevertheless, we think that such the work can only give some insight for " shadows" 
of the real chaotic behavior. 



IV. IS THE KOOPMAN'S CONSTRUCTION THE ONLY ONE POSSIBLE? 

In the previous Section it was argued that in the study of chaotic systems we should 
restrict the class of dynamical systems to that defined on a compact subset of the phase 
space T. This implies that the algebra of observables associated with such the system is the 
set of all complex valued continuous functions defined on Q where Q is a compact subset of 
T. Furthermore, equations of motion for the system - the Hamilton's equations (jl|) should be 
restricted to Q. But one can pose a question: what is the proper choice for Q. If we consider 
a conservative system (such systems have the dynamics given by the Hamilton's equations) 
the most natural choice is the isoenergetic surface. So, in the remainder of this section we 
assume that Q is a compact isoenergetic surface of some dynamical system. Consequently, 
Q will be also considered as a compact topological space. We denote by C(Q) the Banach 
space of all C- valued continuous functions on Q. C(Q) with pointwise defined multiplication, 
complex conjugation and the supremum norm is an abelian unital C*-algebra. By the way, 
we stress that the compactness offl is "translated" in the algebraic language by the property 
of C(Q) to be a unital algebra. But now we have the following problem: C(Q) has the very 
rich algebraic structure, in particular one can write a square. So the concept of nonlinear 
map on C(Q) is perfectly well understood (cf. Example 2, subsection VC). On the contrary, 
the situation on Q is much less clear. To see the problem let us consider, as an example, 
one dimensional harmonic oscillator. Then Q in this case is just the ellipse. So, a natural 
question is emerging: Is it possible, starting from Q, to define such an algebraic structure 
W which is compatible with that of C(Q) ? In the just mentioned example of a harmonic 
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oscillator, the ellipse can be naturally equipped only with topological and measure structures 
and these structures are quite enough for the standard Koopman's construction. However, 
the ellipse has not an algebraic structure and therefore it is impossible to study a nonlinear 
lifting in this setting. To clarify this question let us emphasize that, here, the nonlinear 
lifting is understand as the procedure leading from a nonlinear map defined on Q C T to a 
nonlinear map defined on the set of all observables. To solve the problem we propose to use 
the construction which was called "the nonlinear version of the Koopman's construction" 
(see |[43|| ). To describe the solution of our problem we need some preliminaries. Let us 
recall that a pure state tp on an abelian C*-algebra A (called also a character) is a linear, 
continuous, positive and normalized functional on A such that p(ab) = p(a)tp(b) for a,b G A 
(cf. Appendix). Denote the set of all characters by V. Let us add that by Gelfand-Naimark 
theorem we can identify an abelian C*-algebra A = C(Q) with the algebra of all continuous 
functions on V where V is equipped with the pointwise convergent topology. Moreover, one 
can establish one-to-one correspondence between Q 3 x ^ p £ V. So to keep notations of 
the previous sections we shall identify Q with V. Let Wq denote the algebra of continuous 
functions on C(Q) generated by the set V(= fl) of all characters of C(Q). Elements of 
Wq are of the form $ = W((fi, ...,(p n ), where n G IN, W is a polynomial of n commuting 
variables with complex coefficients, and ipi, ...,(p n G V. Using the correspondence between 



the set f2 and the set of all characters (cf. |43[], see also |3T]| for details of the theory of 
operator algebras) one has the following inclusion 

tt C Wn- (17) 

Now, let us remind that the Gelfand transform is defined as a map 

~:A^C{V) (18) 

which to every a G A assigns a function 

3(v) =¥>(«)■ (19) 

where ip is a character. Consequently, C(V) can be considered as the image of A under the 
Gelfand transform. 

Let a G A. We define the following complex valued function a on Wq: 

a(W(<p!, Pn)) = W{pt{a), . . . , p n {a)). (20) 

It is easy to show that a is in fact an extension of a over Wn. Namely, in the case n — 1 in 
(POD one has 

~a(W(p)) = W(p(a)) = p(W(a)) = W(a)(p). (21) 

Here the second equality follows from multiplicativity of the character p. Further, we shall 
consider dynamical maps on Wp which are of the following form 

9(W(p u ...,<p n )) = (T(H/))(r(^), • • • , r^)), (22) 

where T is a polynomial of one variable, and r : V —>■ V is a continuous map. Let us denote 
the set of all maps in this form by 0. Now we are in position to give (cf. 
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Theorem 1 For an arbitrary map 9 G there exists, in general nonlinear, the well defined 
map 

C(V) 3a^U e ae C{V) (23) 
and 9 is uniquely determined by Uq where 

U e = U T o T. (24) 
U T is the operator defined by (U T f)(x) = f(r(x)). 

Thus, we got the affirmative answer to the question posed in the beginning of this Section. 
There is a possibility to introduce the algebraic structure Wq in such a way that the mul- 
tiplication in this structure is obtained from that of A. Therefore, nonlinearity of a map 
on Wq is compatible with the nonlinearity of the map on A. This gives a possibility of 
nonlinear lifting of a map on Wq to a map on A. 

To make the basic idea of the presented construction more clear let us again consider the 
example of one dimensional harmonic oscillator. Then, the set Q is equal to the ellipse. Wq 
can be considered as "an algebraic reconstruction" of the phase space T (from the singled out 
subset Q) in such a way that its multiplication is compatible with the natural multiplication 
of the algebra of (classical) observables C(Q). The main difficulty in carrying out this 
reconstruction is that there is still one question unanswered. Namely, we do not know a direct 
relation between T(D Q) and Wq. On the other hand, this nonlinear version of Koopman's 
construction suggests rather strikingly the explanation why among the well known examples 
of chaotic systems the periodically perturbed ones play a prominent role. Such the systems, 
in general, are not conservative ones. Let us consider the one dimensional kicked harmonic 
oscillator; a nearly paradigm of chaotic system. Let Q be again the isoenergetic surface of 
the harmonic oscillator. To build Wq we should construct over each point of Q a "fibre" 
which is a free ring of one variable. Then one can suppose that the free motion is determined 
by a map r defined on the isoenergetic surface, i.e., r : Q — > Q. The Theorem [I] says that r 
can be linearly lifted to a map on A. Then a kick (perturbation) acts on the system. Within 
this picture, one may say that a kick corresponds to a motion along a fiber. Consequently, 
a kick is reflected on the set of observables as a nonlinear map (see Theorem |1|). Thus, we 
have got the explanation saying that kicked systems can be related with the "algebraic" 
nonlinearity. We remind, a nonlinearity is the necessary condition for a chaotic behavior. 

Having clarified the basic concepts for classical dynamical system we want to pass to the 
basic subject of the paper, namely to the problem of existence and description of quantum 
chaos. 

V. QUANTUM MECHANICS 

It was pointed out in Section II that the basic concepts for a description of a classical 
dynamical system are: algebra of observables, states, and prescription for the equation(s) 
of motion which can be formulated in one of two pictures (equivalent in the sense of equal- 
ity (|9])). To quantize the (classical) scheme it is enough to change the realization of the 
algebra of observables with the relevant modifications of description of state and evolution. 
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Therefore, we replace the algebra C(Q) (where Q is a compact subset of T) by a unital (non- 
abelian) C*-algebra A. We repeat: the unitality of the algebra is the algebraic translation 
of the topological compactness property! Then, the states are linear positive, normalized (so 
continuous) functionals on A. The time evolution is described by quantum counterparts of 
the Liouville and Hamilton's pictures, i.e. by the Schodinger and the Heisenberg picture 
respectively. Our objective in this Section is to convince the reader that the Heisenberg 
picture is the suitable framework for investigations of quantum chaotic problems. 



A. The Schrodinger picture 

1. Linear Schrodinger equation 

This is the quantum counterpart of the Liouville picture for a pure state. The basic 
equation of motion is 

. h d^(t) .^d^(t) . . 

i — = ih — — = H^(t) (25) 

2tt dt alt W y 1 

where \l/(t) denotes the wave function. Its solutions are of the form 

ty(t) = e-^^(0) (26) 

Thus 

= ||^(0) -tf'(0)| | (27) 



We observe that the equality ( |27| ) implies that there is not any growth of separation of 
orbits which at time zero differ by a small vector! Thus there does not exist a room for 
a sensitivity on initial conditions. Consequently, there is no chance for positive Lyapunov 
exponents. We want to emphasize that exactly the same situation we have for a classical 
dynamical system with dynamics given in the Liouville picture. Not taking into account the 
last remark, the long debate on existence of quantum chaos has led to a conjecture that 
quantum mechanics suppresses chaotic behavior in such a way that there is only a room for 
some "shadows" of chaos. Adopting this point of view one could only study special features 
of a wavefunction which describes quantum counterpart of classically chaotic system and to 
investigate some properties of the correspondence between classically chaotic system and its 
quantum counterpart. In other words we have the quantum analogy of the search of chaotic 
signatures within the Liouville picture (cf. the last paragraph of Section III). 

Let us briefly discuss a typical result in this direction. The most famous proposition 
concerning a characterization of special properties of quantized classically chaotic systems 
is given by the following conjecture: There is a correspondence between classical (regular 
or irregular) dynamical system and the energy level statistics of the Hamiltonian of the 
quantum counterpart of the considered system. Namely, 

regular (classical) system <-> Poisson distribution (28) 
irregular (classical) system <-» Wigner — Dyson distribution (29) 
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where the Poisson distribution is given by P(AE) = ae bAE : a, b are some positive constants 
while the Wigner-Dyson distribution is given by P(AE) = A\AE\ a e-^ AE ^ 2 , with A,a,/3 
some positive constants. However, recently (see |]18|) it was shown that there are infinitely 
many counterexamples to the above conjecture. Thus, although there are some computer 
experiments verifying this conjecture for some models it is difficult to say that there is a 
well defined correspondence between energy level statistic and the nature of evolution. 

As a matter of fact the same picture is emerging from the analysis based on the other 
apparently basic features of the wave function of quantized (classically) chaotic system. To 
specify them let us list: scars, revivals, localization of wave function, level repulsion leading 
to regions of complicated avoiding crossing. To the best of our knowledge, there does not 
exist any rigorous proof of a correspondence between such features and quantum chaotic 
systems. Therefore, looking for an answer to the question whether there is "true" chaos in 
quantum mechanics we drop this way of investigation of the posed problem. 



2. Non-linear Schrddinger equation 
Let us consider the Schrddinger type equation (in general, nonlinear) 

1%^- = H i; + V{if>) (30) 
at 

where ip G TC, Ti is a separable Hilbert space, Hq is the free Hamiltonian (yielding a linear 
evolution), and V(ifj) is a (nonlinear ) potential which is defined below. In other words, we 
begin with a field-theoretic model where the states at a given time are represented by vectors 
in TC while the nonlinearity enters through their time-dependence, i.e. the nonlinearity is 
described by (|30|) . To give the definition of V(-) let us denote by {e^} a basis (independent 
of t) in TC. Define 

H3 ^ ^U(iP) = {< e k ,iP >} 6 4 (31) 

{ti is the Hilbert space of sequences {cti}, ctj a complex number, such that J2i |««| 2 < oo.) 
Denote by U the unitary map from TC onto li. We assume that a function v({oti}) = f {v(ai)} 
satisfies the following condition 

\v («j)| 2 < oo provided that ^|ctj| 2 <oo (32) 

i i 

The nonlinear term V(ip) in flSTf ) is defined by 

Vty) = U~ l v(Um (33) 

Let us remark that if one takes another basis, say {ft}, in TC then (|3l|), (|32f) , and ( |33l) give, 
in general, another nonlinear operator V'{ip) on Ti. However, one has 

W~ l V'{Wil)) = V(ip) (34) 

where W is the unique unitary operator on Ti such that ^({6^}) = {/&}■ Thus the operator 
V(ip) is defined in a "covariant" way. 
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As examples of nonlinear potentials leading to nonlinear type Schrodinger evolution (p0|) 
we can mention V\(i)j) = ifjipifj (its Heisenberg counterpart has applications to molecular 
dynamics and nonlinear optics, |19[ see also P5[ , or V^VO — ^(VO f° r V* C 7i 



where the domain H of the (unbounded) operation log(-) is determined by the condition 
(such evolution describes frictional effects in dissipative systems, [j3"5] , Other 



examples of nonlinear Schrodinger equations can be found in ( |24| , [ |13|1 , [34|) 



However, one may ask whether the equation is of the fundamental nature or it is a 
kind of approximation. Although there are some arguments that nonlinear Schrodinger type 
equations could be fundamental ones (see e.g. see also the argument given in Subsection 



VB), this question is still unanswered. Namely, one can also consider this type of equation 
as a result of an assumption that a part of a larger quantum system was regarded as "heat 
bath" . The lack of convincing arguments can be also related to the question of the nature of 
frictional effects. If one admits the conjecture that the Newton's dynamics is more general 
one than that of the Hamilton's the answer to the just raised question could be affirmative, 
i.e. one can expect to find "genuine" nonlinear quantum maps. Let us add that, in such the 
case, even on the level of classical mechanics the Liouville picture based on the Hamilton's 
dynamics would be not the most general prescription for the evolution in the sense that 
nonlinear corrections would be necessary ingredient of the equations of motion. 

Nevertheless, having a nonlinear Schrodinger equation one can ask for a possibility of 
chaotic behavior. To answer this question one can proceed as follows. We begin with 
rewriting (|30| ) in terms of Fourier coefficients, i.e. using the basis {e k = \k >} in 7i and the 
formula (0) we arrive to 

ih— < *# >\k>=Y,< k\H i/j >\k>+Y,< k\V(il>) > \k > . (35) 

k k k 

Now, just for simplicity, let us put H = 0. Moreover, to give as simple as possible arguments 
we take a concrete form of the nonlinear potential. Thus, in the definition of V(-) let us 
put v(a) = KaJ2n${t ~ n ^o) + 2® aa wnere K an d X are positive constants, to is a fixed 
period of time. Such potential was designed to describe the nonlinear evolution in the Kerr 
medium with kicks where the kicks are defined by the first term in the formula for v, i.e. 



KaJ2n$(t ~ n to) (for details see []45|| , |J46| ). Hence 



ih—z k = Kz^^2S(t-nt ) + -z^z k z k k = 1,2, ... (36) 



where z k =< k\ip >. Let us note that owing to the properties of Dirac's S, fl3"6j ) can be 
considered as a system of two equations; the first one determines the dynamics in the period 
between kicks while the second one gives the parametric process describing kicks in the time 
t n = nt , n — 1)2, ... (cf. [^6|). Clearly, ( p6|) is the (standard) nonlinear differential equation 
and an application of classical methods leads to solutions of (^) which admit a positive 
Lyapunov exponent (cf. |46|]). In other words ( |30D leads to an example of nonlinear 



Schrodinger type equations with sensitive solutions with respect to initial conditions where 
the Lyapunov exponent, A^, is defined as A^ = lim^oo jlog^^^ and is expressed in 
the Fourier coefficient terms (cf. [0). So we conclude that within the nonlinear type of 
Schrodinger evolution there is a room for a chaotic behavior. 
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3. Density matrix formalism 



The genuine quantum counterpart of the Liouville picture is that given in terms of density 
matrices. So, let us consider the quantum Liouville equation, called the von Neumann 
equation, which is of the form 

q = ~[H,q] (37) 

where H stands for the Hamiltonian of the considered system, g is a density matrix while 
[■, ■] denotes the commutator. Clearly ( p7|) is a linear differential equation and one can ask 
for nonlinear version of this equation. We recall that nonlinear equations for time evolution 
of density matrix are known and used in physics (cf. ||, ||10|| , @). Furthermore, similarly as 



for the Schrodinger equation, the question concerning the nature of the equation (|3T|) as well 
as its nonlinear corrections can be posed. Basically one has the same answers. Therefore 
we will not discuss this point here; let us pass to an illustrative example instead. Let us 
consider a dynamical system with dynamics given by Hartree - type evolution equation (cf. 



^ = -i[*to,*] (38) 



Pt\t=0 — Poi 



(39) 



where H(p) = [Tr(pQ)]Q with Q = Q*, and g is a density matrix on a finite dimensional 
Hilbert space TC. We recall that this type of equation can be obtained in the mean field 
limit for interacting quantum systems. The explicit solution of ( j38f) is given by 



$t(Po) = exp [-^-Tr(Qp )Qt]p exp [^Tr(Qp )Qt). 



(40) 



Clearly, the solution $(•) of (|38D is a nonlinear quantum map. 

To summarize, within the quantum Liouville picture there are examples of nonlinear 
evolution equations. Consequently, this can be taken as a motivations for our investigations 
of stability of dynamics determined by such the equations. 



B. The Heisenberg picture 



We recall that this is the quantum counterpart of the Hamilton's picture (the picture 
from the classical mechanics, cf. Section II). The fundamental equation of motion is of the 
form 

where A is an operator in the algebra of observables, while H stands for the Hamiltonian of 
a system. The solution of fl4"Ip can be written in the form 
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A h-> e% Ht Ae~% m 



(42) 



Our first remark is that the Heisenberg formula for equations of motion is in harmonious 
relation with the very rich algebraic C* structure of the matrix formulation of quantum 
mechanics. On the contrary, the Schrodinger picture does not have this property: a Hilbert 
space has only a linear space structure while the set of density matrices for a real physical 
system is a Banach *-algebra. There is only one exception. Namely, for finite dimensional 
models of quantum mechanics, the set of trace class operators is equal to the set of all 
bounded operators. Therefore, for such particular models one can also use the C* algebraic 
technique for an analysis of evolution of density matrices. The next remark is related to a 
possibility of getting nonlinear operator equations. To be more specific, let H be a Hamil- 
tonian of a physical system such that H = H(A, B) is a function of several noncommuting 
dynamical variables. Obviously, in general, the definition of a function of two (or more) 
noncommuting variables is not a clear stuff. However, in Quantum Mechanics there are 
many examples of Hamiltonians of this type which have arisen as a result of the quantiza- 
tion procedure. So, here, we can put for the concrete H(A, B) the relevant Ansatz. Then, 
the solution of the Heisenberg equation of motion for A (the equation needs not be linear in 
A) 

A ^ e iH ( A,B)t Ae -fH{A,B)t (43) 

does not have to be the linear function (in A). Clearly, we have not such a possibility for 
the Schrodinger evolution of pure states. 

Let us consider the question of nonlinear dynamical maps in the Heisenberg picture 
in detail. As it was mentioned at the beginning of this Section the Heisenberg picture is 
quantum counterpart of the Hamilton's picture. Let us remind (cf. Section II) that the 
Hamilton's equations also lead to a one parameter group of linear maps on the algebra A 
of observables where A = C(Q). In particular, the Poisson bracket is nothing else that the 
derivation on A (cf. Appendix). However, as it was pointed out in remarks following (^) 
there is a possibility, for the very special choice of observables, to study of time evolution 
without the linear lifting. Namely, the choice of {q} and {p} can lead to nonlinear differential 
equations with solutions given by nonlinear maps (even to the solutions which are sensitive 
with respect to initial conditions). Now let us turn to quantum mechanics. We know, e.g. 
from quantum optics, that the creation a) and the annihilation a operators can be treated 
as a substitution for q and p. Let us take this point of view and let us put in ([EJ, H3) 



A = a' , B = a. Then, we can expect to obtain nonlinear dynamical maps for some choice 
of Hamiltonians. Consequently, we can expect (similarly as in classical mechanics) to find in 
this way nonlinear equations and maps with interesting (non)stability properties. We want 
to add that another way of introducing nonlinear maps in quantum physics will be studied 
in the next subsection. 

We close this subsection with a brief discussion on the equivalence of Schrodinger and 
Heisenberg pictures. Similarly, as for the classical mechanics (see (^)) we have 

< A > av (t) =< A > m = (tt(t), Att(t)) = (e-i H{A ' B)t y(0),Ae-i H{A ' B)t V(0)) (44) 
= (¥(0),e*"< A '^4e-* ff <^¥(0)) =< A{t) >* (0 )=< A(t) > av 
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To get a better understanding of the nature of nonlinearities described by (f43|) let us replace 
A in 



A h+ e iH(A,B)t Ae -±H{A,B)t (45) 

by < A > av . This can be done, e.g. by performing some approximations in averaging 
procedures (cf. [pTjl , |45| ). In this way we can arrive to 

< A > av ^ /(< A > av ,t) (46) 

where /(•) is, in general, a nonlinear function. Thus, we can obtain the "classical" nonlinear 
maps and we can expect to find examples having sensitiveness with respect to initial condi- 
tions. In fact such examples were found, e.g. see ( pTfl , HJ) Obviously, in our discussion of 
stability properties of quantum systems we want to avoid the obscure procedure leading to 
([IB]). To this end we shall introduce quantum Lyapunov exponents (cf. Section VI). 

To finish our discussion on equivalence of both pictures we should say few words on the 
relation between the nonlinear Schrodinger equation ([3(]) and the Heisenberg equation ( f41"| ) . 
Obviously, the prescription (|44] ) based on bilinear form (•, •) (given by the scalar product) is 



not the relevant one. However, to pass from a differential equation for a state to an operator 
equation involving an observable we can use "the second quantization rule" (cf. ||14|| , [ |I~9[ 



sec 



also p9|| ). To explain this idea let us restrict ourself to the model described by equations 



(|35l) and ([36]). Then, the equation of motion of the annihilation operator (respectively 
the creation operator a*) follows from the equation of motion ( |36"|) in which we replace the 
Fourier coefficients and HI (so complex numbers) by the boson creation and annihilation 
operators a { and a* (so operators). Thus, we have 



^Zu a i = Ka *j ^ ~~ n ^°) ^ — a *j a i a i (47) 



■j t a i = Ka j s ( f - nt o) + \ 

The remark following (^fj) is also applicable to this case. Let us observe that ([47D follows 
from the Heisenberg equation of motion 

with H = z| Y,j[( a j) 2 ~ ( a j) 2 } J2 n d(t ~ n to) + f J2j a *j a ] a j a j- Thus, the second quantization 
methods links nonlinear Schrodinger type equation with the undoubtly fundamental equation 
of motion of the Heisenberg form. Moreover, the latter can lead to nonlinear operator-valued 
equations. 

To summarize this subsection, we conclude that within the Heisenberg picture there is a 
room for quantum counterparts of "nonlinear" type evolution. 



C. Are there many nonlinear quantum maps? 

We have pointed in previous Sections that although nonlinear relations are not sufficient 
for chaos of classical systems, some form of nonlinearity is necessary for chaotic dynamics. 
We observed that the Heisenberg equations of motion as well as the (nonlinear) Shrodinger 
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equations allow some form of nonlinearities. Besides we would like to know whether a non- 
linearity can arise from the "nonlinear" quantization of nonlinear systems. We remind that 
in the theory of classical chaos we frequently study the dynamical maps alone. This arises 
the question for a quantization such the dynamical maps without referring to Hamiltonian 
equations. In other words we are talking about a possibility of quantization of Newtonian 
systems. Working in this direction, we have shown (cf. Section IV) that there exists the 
nonlinear version of Koopman's construction. Here, we indicate how this construction can 
be used to argue that the well known class of nonlinear completely positive maps defined on 
a C*-algebra (see Section VI and the Appendix) may be of the "physical" nature; for another 
motivation to study nonlinear completely positive maps in quantum dynamical systems see 

Let again A a bei = C(Q) = C(V) be an abelian C*-algebra with unit where Q is a compact 
(Hausdorff) space (cf. Section IV). Within the algebraic structure Wn we can consider the 
map ip I— >■ (p n . We observe 

a{ V ) = ~a( V )^~a( V n ) = {^)( V ) (49) 

where a G Aabei- In other words, there exists the correspondence between the map (p i— > tp n , 
<p G Q, (p n G Wfo and the map a i— > a n , a, a n G Aabei- Consequently, if a subset Q of the 
phase space is algebraized properly, i.e. it is endowed with the Wq structure we are able to 
translate nonlinearities of maps defined on Wfo to that defined on A a bei- Let us illustrate 
this point. 

Example 2 Let us define 

C d f = rf(l-f) (50) 

where 1 stands for the identity function, /(■) G C(V) while r is a constant. We have 

C d M = (rf- rp)( V ) = Mf) - Mf) = rcp(f) - r<p 2 (f) = (51) 
- (r(p - <P 2 )) (f) = (Cp){f) (52) 

where tp G T 3 , / G C(V) and finally f is the element of A which defines the Gelfand transform 
f . Consequently, our example demonstrates rather strikingly that the logistic map (one of 
the most famous "chaotic" maps) can be considered both in Wp and C(V) terms. Evidently, 
C d as well as £ are nonlinear maps. 

Let us pass to non-abelian case, i.e. we replace A a bei by A, where A is an arbitrary 
nonabelian C* algebra. We remind that A can be considered as a noncommutative analogue 
of the space of bounded continuous functions, i.e. bounded continuous functions over a 
"quantum plane". Here, we are using the phrase "quantum plane" in the sense of Manin 
(see |Q). Let us repeat the argument given at the beginning of this subsection but now 
applied for the quantum plane. We can expect that the maps of the form 

a i— > c n a n (53) 



15 



where a G A, c n is a number, have arisen as a result of lifting nonlinear maps defined over 
the quantum plane (cf. Theorem 1, Section IV). In other words, it is our hypothesis that the 
maps of the form $3Q, or more generally completely positive nonlinear maps on a C* -algebra 
can be considered as a result of (nonlinear) lifting nonlinear maps defined on the quantum 
plane. This would mean that the quantization of a Newtonian system could lead to genuine 
nonlinear quantum map. 



VI. QUANTUM LYAPUNOV EXPONENTS. 

In this Section we develop the theory of quantum Lyapunov exponents (cf. f4"0|| , [fflj). 
Let A be a C*-algebra with the unit 1, r : A — > A a quantum map, i.e. 

1. r is positive: r(A*A) > for all A G A, 

2. r(l) = 1. 

The pair (A, r) consisting of a C*-algebra A and a quantum map r will be called the quantum 
system. We propose the following generalization of the classical Lyapunov exponent: 

Definition 3 Let (A, r) be a quantum system where r is differentiate (in the Frechet sense) 
map. Let D x r n (y) denote the directional derivatives of r composed with itself n times at a 
point x in a direction y. Then, the limit 

X^r;x,y)=\im ^log\\(D x T n )(y)\\ (= A 9 ), (54) 

whenever exists, will be called the quantum exponent. 

Let us comment on this definition. Firstly, it is clear that ( |54| ) can be generalized for 
continuous dynamical systems in an obvious way (see Section III). Secondly, one can study 
the following weaker version of quantum characteristic exponent: 

Xl(r;x,y)=lim Uog\»{{D x T n ){y))\ (= A 9 ), (55) 

where /i is a state on A. However, as we are basically interested in quantum chaos it is 
enough for us to look for the largest Lyapunov exponent. As X q plays such the role (see 
properties of quantum exponents) we shall restrict ourself to definition involving the norm. 



Thirdly, the noncommutative generalization ([54]) of characteristic exponent is not the only 



one possibility. Namely, within the Connes's noncommutative geometry, ||21|| , the required 
differential structure can be introduced by derivations (see Appendix for definitions) {5} of 
the algebra A associated with a physical system. Then using basically the same idea (see 
f26fl ) one can define 

A" = lim -log||(5n(x)|| (56) 

t >oo £ 

where x G T>(5) and r t stands for (continuous) dynamics of the system. Moreover, we should 
mention that the presented two examples of definitions of quantum Lyapunov exponents do 
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not exhausted all possibilities (cf. fl53| , |56| , |57|). However, studying the properties of 



quantum exponents as well as looking for quantum counterparts of chaotic systems (see 
Section VII) we have found X q as the best suited. Finally, the slight modification of the 
definition of Lyapunov exponents described in Remark 1 (see Section III) can be also done 
for quantum exponents. In fact, we shall use such the modification in our examples (5) and 
(6) (see Section VII). 

Let us turn to the question of existence of the limit in the definition of X q . In other 
words, to present a convincing argument that X q is well defined we should give an example 
of sufficient conditions implying the existence of the limit in (0). Let us assume: 

1. r is a completely, in general nonlinear, positive map (see Appendix for appropriate 
definitions). In fact, this is the most important assumption for the next theorem. 
Namely, this condition implies a smoothness of the quantum map r; even one can 
say that the complete positiveness implies a kind of Laurent expansion for r, (see 
Appendix for details). 

2- ||t'(0)|| < C\ for all I G N and some positive Cj, and 

6 T = {x + : \r\x)\ < C 2 ||x|| + 11^(0)11} ^ (57) 

for some positive C*2 and all / G N. 
3. Finally 

\\D x r k (y)\\ >C k (x,y) (58) 

for some positive C(x,y) and all large k G N. 
Under the above assumptions one can prove (cf. |4lH ): 



Theorem 2 Let r : A — > A be a map such that the assumptions (1), (2), (3) are satisfied. 
The following limit 



lim - log \\D x r n (y)\\ (59) 

n^oo 77, 



exists for x G T . 



The just given theorem clearly shows that the notion of quantum characteristic ex- 
ponent X q is well defined for a large class of dynamical systems. Therefore, having the 
well-established notion of quantum counterpart of Lyapunov exponents X q let us turn to the 
description its basic properties. We have proved (see P2|) 



X q (x, y) = X q (x, ay) for aGR\{0}. 
Note that since the map y —>■ D x r n (y) is linear one , it is natural to set 

X q (x,0) = -oo. (60) 
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Let X q (x,y) > X q (x,z) > — oo and, additionally, let r satisfy assumption (iii) of the 
theorem in the direction y + z. Then: 



X q (x,y + az) < X q (x,y), (61) 



where a G R. 



The function y — > X q (x,y) as the limit of continuous functions (in y) is, in general, 
the Baire function of type I. In particular, the set {y\y — > X q (x,y) is a continuous 
function} is dense in A. 

Let A be the C*-algebra generated by a fixed self-adjoint operator and identity on 
some Hilbert space, i.e. A is an abelian C*-algebra. Furthermore, let r : A — » A be a 
smooth map while a state on .A. Note that such dynamical system (A, r, 0) should 
be considered as the classical one. Then, if some mild technical conditions are met (cf. 
[12]) Definition^ leads to the largest classical characteristic exponent. 



It is important to observe that the above listed properties of X q are reminiscent to that 
of classical characteristic exponents (cf. ||49|| ). Therefore, we conclude that the notion X q 
is the well defined quantum counterpart of (the largest) characteristic exponent. 



VII. EXAMPLES. 

In this Section we want to study concrete models of dynamical systems. In particular, 
we are interested in examples of systems with irregular time evolution, i.e. with systems 
with positive quantum exponent. Moreover, we shall consider examples with analytically 
calculated exponents. Unless otherwise stated the existence of quantum Lyapunov exponent 
for each example was proved separately. Let us begin with the simplest model. 

(1) A dynamical system composed of A-level quantum system with dynamics given by 
a linear dynamical semigroup St of contractions. The assumption of local nonexpansiveness 
is always satisfied but the condition of variability is satisfied for some directions. In general, 
X q (x,y) G (—oo,0). As an illustration let us consider S t x = e~ M x, where t > and A > 0. 
Then X q (x,5x) = —A. This means the stability of the considered semigroup dynamics. 

(2) A dynamical system with dynamics given by Hartree - type evolution equation (cf. 

0): 

±p t = -i[H(pt),pt] (62) 



pt\t=o = po, (63) 

where H(p) = [Tr(pQ)]Q with Q = Q*, and g is a density matrix on a finite dimensional 
Hilbert space 7i. Let us recall that this type of equation can be obtained in the mean field 
limit for interacting quantum systems. The explicit solution of ( |S2]) is given by 

$ t (p ) = exp [-iTr(Q Po )Qt] Po exp [iTr(Q Po )Qt}. (64) 
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Thus taking the discrete time in fl64|) and putting r = $1 in (|54D we get: 



lim — log || — i[Tr(5pQ)Qn, p] + 8p || = 0. (65) 



In (65) [•, •] denotes the commutator. Clearly, assumptions of our theorem are satisfied. As 
X q = 0, it follows that Hartree-type evolution ( |62D is an example of the regular motion. 

(3) As the next example let us consider the dynamical system (£>(H),$) where B(H) 
denotes all linear operators on a finite dimensional Hilbert space while the map $ is given 
by 

$(p) = p 2 (66) 

for p* = p from unit ball in B{7i). We want to stress that this example is rather "pure" 
mathematical one in the sense that it is well defined mapping as well as the "second order" 
term in the expansion of completely positive nonlinear maps on C*-algebra (cf. Appendix). 
But, it is difficult to provide a (physical) Hamiltonian leading to the considered map. On 
the other hand, let us recall that the "square" is the basic ingredient of the logistic map - 
the famous paradigm of (classical) chaotic systems. So, if one assumes that the nonlinear 
version of Koopman's construction can be applied to nonlinear quantization (cf. Section IV 
and subsection Vc) one can expect to find the physical interpretation of this map. Finally 
we want to add that the transformation of this type was used by Gisin and his coworkes 
(cf. ||) for the study of the theory of quantum computers. So one can conclude that such 
systems may be of great significance for a general theory. 

Having the clear motivation let us turn to an examination of its properties. It is easy to 
observe (cf. Euler theorem) that 

||D p $"(p)||=2"||p|| 2 "- 1 (67) 

for p 7^ 1 and p ^ Projectors(7i). To prove (p?[), it is enough to note that \\p\\ = 
sup\£ Spec trum(p)\M- Clearly, the assumptions of the theorem are satisfied. Consequently, 
(|67|) leads to: 

A"(p,p)=lim - ((2™ - 1) log || || + n log2). (68) 

As a conclusion, this mathematical example shows that the quantum exponent X q can also 
be positive (put p : \\p\\ = 1, then X q (p,p) = logl > 0). Thus the quadratic map, in some 
regions, can exhibits the irregular form of the "motion" . 



In the next three examples we reexamine dynamical systems which were semiclassically 
treated recently. Namely, using a kind of approximation (cf. our discussion following ([EJ), it 
was shown that these models of quantum origin display signatures of chaotic motion. Now 
we want to treat these systems in pure quantum way and compare results. We want to 
add that in these examples the Hilbert space Ti of a system is not assumed to be of finite 
dimension. 
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(4) Let us consider a quantum two-level system interacting with a single mode of the 
electromagnetic field (cf. p7 |, ||32|| ). The Hamiltonian for a such system can be taken as 



H = -hu) a z + tiu(ata + -) + hXoa x (a + at) (69) 

where a x , a z are Pauli matrices and and a stand for bose creation and annihilation 
operators. The energy separation of the two-levels is equal to Tiloq while the frequency of 
the radiation mode is u. The time evolution of the system is given by the following operator 
equation 

A = j:[H t A] (70) 

for an observable A, where the dot denotes a time derivative. Let us remark that the 
Heisenberg equation fl7D| ) gives for this model first order coupled nonlinear operator equations 

(cf. my- 

a x = -uay (71) 

& y = oo a x — 2Xa z (a + at) 

a z = 2Xa y (a + a') 
(a + a^) = —iuj[a — a') 
(a — at) = — iu{a + at) — 2i\a x 



Although the Hamiltonian H depends, in general, on the time t, one can solve (|70|) and 
subsequently compute the quantum exponents X q (see f32|| ). In particular one has: 

X q (a u a J ) = (72) 

for i,j = 1, 2, where 0\ = a x and o-i = a z . Consequently, this quantum optic model exhibits 
quantum regular evolution evolution. Let us add that using some approximations a signature 
of chaos was shown (cf. [|27j1 ). 

(5) This example is the quantum counterpart of a parametrically kicked nonlinear oscil- 
lator (cf. [j33|). Its hamiltonian H is of the form 

H = ^(a t )V + ih-[{at) 2 -a 2 }- £ 5{t - nt Q ) (73) 

^ n=— oo 

where \ and n are constants characterizing the system, at and a stand for bose creation and 
annihilation operators, and to is the period of free evolution (i.e. the evolution described 
by Hnl = f (a^) 2 a 2 ). Let us remark that H^l gives the nonlinear part of the evolution. 
Obviously, S(-) stands now for the Dirac distribution. Then, the Heisenberg equation of 
motion, 

A = UH,A] (74) 



20 



can be solved in terms of self-adjoint operators $ and II, where $ and II are defined by: 
a = $ + iIT. The solution of the Heisenberg equation for the considered model can be written 
as: 

(*(t+)\ =m ( e^cos^o e'-sin/iSo \ f$(C-i)\ (7 ^ 
\U{t+)J \-e- r smfiB e' r cos fiB J \U{t+^) J 1 ; 

where n = xto, r the effective constant for kicks, and Bq = a^a — \ = ($ 2 + II 2 — 1). 



Let us remark that ( |75j) gives a nonlinear evolution since the matrix on the right hand 
side of ([75|) is the nonlinear function of operators $ and U. We want to study the stability 
properties of the evolution ([75]) . We restrict ourself to an examination of the Lyapunov 
instabilities of the quadrature components of the electric field during the time evolution 
given by the formula (|75|). To do this let us define the quadrature operators: 

<F = -[e ie a + e- ie at] (76) 

and 

n e = -[e ie a-e-V] (77) 

where e G [0,27r]. Let us note that the operators $ e and IT are related to the amplitude 
components of electric field. At this point in order to avoid any confusion some explanation is 



necessary. Namely, the Hamiltonian (73) leads to the nonlinear evolution map ( |75|) . However, 
from physical point of view, an examination of the stability properties of dynamical map does 
not allow any change of annihilation (so also creation) operators nor any modification of 
the fixed Hamiltonian. This makes the big difference between the treatment of classical and 
quantum models. Namely, in classical models to study stability properties of dynamical 
maps on the phase space (for the special choice of observables, cf. Sections II and III) we 
can take a slightly different coordinate and momentum. Here, creation and annihilation 
operators play the role of coordinate and momentum. Again, we made a very particular 
choice of observables. The result was nonlinear operator-valued equations. Therefore, to 
find physical examples with positive quantum exponents X q and to keep the uniqueness of 
creation and annihilation operators we are examining the quadrature operators. In other 
words, the quadrature operators are examples of dynamical variables in our model which 
can be varied without any change of a (a') and the Hamiltonian. Thus we shall study 
quantum exponents in the sense of Remark 1 (see Section III). 

To compute quantum exponents for quadrature operators in the considered model we 
should find an explicit formula for the norm of ||D e r n (y)||, where in the example considered, 
by r we denote the time evolution of the system between two successive kicks as well as 
the effect of the first kick. The variable y stands for an element of the set {$, n}. Clearly, 
r n denotes r composed with itself n times. D e is the derivative with respect to the phase 
angle e which describes a "rotation" in ($, n) variables and which is used in the definition 
of quadrature operators. But, since a* and a are unbounded operators, in order to com- 



pute ||D e T n (-)|| we should introduce ($,5,^), the cut-off in the ($,n) variables (see |33 



for details). Then, replacing the original dynamics by its well defined approximation and 
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changing the variables ($, II) to ($,11) one can find such set of parameters (x, k, to) that 
A 9 (ft) > 0. It is remarkable that this result does not depend on the performed cut-off. In 
other words, for large enough r (r is the effective constant for kicks) and some values of 
[i = ton, the quantum characteristic exponent X q for quantum variable 11,5 is strictly posi- 
tive. To comment on this result let us note that in polyparametric cases physics as well as 
mathematics allow numerous combinations of stability in certain directions and irregularity 
in others. Consequently, in such the cases one can expect chaotic features of the evolution of 
the model for some values of x, r, r, k and regularity for others. Here, we got a confirmation 
of such the behavior (\ q (§s) can be negative for some values of parameters k, x, and r and 
positive for other values of parameters). Moreover, let us add that the approach based on 
some approximations (cf. remark preceding example (4)) to the considered problem gives 



the similar result (see [j4"5] , pESfl ). We conclude that this model exhibits the hyperbolic type 
of instabilities. 

(6) In the last example we consider the squeezed light in a nonlinear medium having 
the second-order susceptibility (x^)- An analysis of such system leads to the following 



equations (see pOJ, section 11): 

j^a(z) = kat(z), (78) 

^-a\z) = ka(z), (79) 
dz 

where a(z) is the annihilation operator, k is the coupling constant and k stands for the 
complex conjugation of k. The dependence of a on z is a result of the one dimensional 
propagation of the electric field along the z-axis. Obviously, the equations (|78|) , (|79|) lead to 

d 2 

— a{z) = \k\ 2 a{z). (80) 

It is easy to observe that the same equation can be derived from the Hamiltonian 

H = hua^a + ih^ 2 -a 2 ) (81) 

where u = 0. Such a Hamiltonian would appear if we take into account an oscillatory 
time dependence of a, and then eliminate the free evolution by the interaction picture. The 
solution to equations (0), ([T9|) is 

k 

a(z) = cosh(\k\z)a(0) + — sinh^l^a^O). (82) 

It is clear that we can apply our definition of Lyapunov exponent X q to analysis of stability 
of the propagating squeezed light. The quadrature operators are taken in the form 

P a W = \ (e ia a(z) + e~ ia a){zj) (83) 
Q a (z) = ^(e ta a(z)-e- ta a\z)) (84) 



22 



where a G [0,27r]. In particular, we have 

P a ( z ) = i([e |fc|2 sina + e- |fc|2 cosa]a(0) 



(85) 



+ 



An easy calculation leads to 



-e |fe|2 sina + e |fc ' 2 cos a] a) 



X q = lim lim - In 

n— >oo z— >oo 2, 



dQ n a (z) 



da 



\k\ > 



(86) 



(87) 



where we put for simplicity 



, fc , — —1, indicates the cut off in Q a (since a and are 
unbounded operators we again have to apply cut-off). 



It is remarkable, see (|87|) , that this result does not depend on the cut-off. Hence, this 
is another quantum optic model with very unstable dynamics. It is important property of 
this model that the classical Lyapunov exponent calculated for the trajectory 



t i ^ (w\Q a (t)\w) 

has the same value as its quantum counterpart; in (| 
a(0)\w) = w\w). 



(88) 



w denotes a standard coherent state 



VIII. CONCLUSIONS. 



We have shown that Quantum Mechanics allows some form of nonlinear equations for 
dynamical maps. Consequently, we have found a kind of coexistence of regular and irregular 
quantum motion; the coexistence resembling very much that of classical mechanics. Both 
forms of nonlinearity, in Classical and Quantum Mechanics, arose as a result of the special 
choice of observables ({p, q} for classical systems, {a, a^} for quantum systems). However we 
want to repeat that there is very important difference in the study of stability properties of 
the corresponding dynamical systems. Namely, studying a stability of classical maps for such 
the choice of dynamical variables one can consider an arbitrary (small) variation of p and 
q. Basically we can repeat this argument for quantum maps defined on the set of states (cf. 
examples 2 and 3, Section VII). However, we can not do this for a and a' ! The proper way 
out of this situation is to study another class of observables, e.g. the quadrature operators 
(cf. examples 5 and 6, Section VII). Namely, although the quadrature operators are built 
directly from a and a\ additionally, they depend on a parameter (phase). Consequently, 
their evolution is described by nonlinear maps while their variation can be implemented by 
that of the parameter (phase). Clearly, this procedure does not effect the uniqueness of 
a, a'. To summarize we have shown that the Quantum Mechanics admits the irregularity 
of the evolution of some variables. We consider the positivity of X q - so irregular motion 
- as the basic signature of unstable dynamics. However, we do not identify this signature 
of chaos with the very chaos of dynamical system since to get a complete description of 
chaotic behavior of dynamical maps some additional geometrical assumptions on the space 
of dynamical variables as well as ergodic assumptions are necessary (cf. Section III). 
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The paper also deals with the related question of derivation of nonlinear quantum dy- 
namical maps without resorting to quantization of equations of motion. We have shown the 
possibility for nonlinear lifting of dynamical maps. These investigations were intended as an 
attempt to motivate our interest in completely positive nonlinear maps on C*-algebras. Our 
result suggests that there is a chance to obtain such the maps via "physical considerations" . 
However this question is at present far from been solved as we do not know the relation 
between the algebraic structure Wo and the phase space T. 

We want to end this paper with the remark that the presented results are steps toward a 
quantum ergodic theory, where a rigorous definition of quantum chaos may be established. 



IX. APPENDIX. 



The theory of C*-algebras can be find in the books of [0, [[5T]] or ||58 . Here we 
recall only few basic results which have been used in the paper. 

• A Banach *-algebra A is called a C*-algebra if it satisfies ||a*a|| = ||a|| 2 for a £ A. 

• Let A be an abelian C*-algebra, i.e. ah = ba for any a, b G A. A character (p, of A, is 
a nonzero linear map, A3 f{a) G (D such that 



Lp(ab) = Lp(a)tp(b) 



(89) 



foe all a,b G A. The spectrum V(A) =V,oiA is to be the set of all characters on A. 

• If (p is a character of an abelian C*-algebra A then f(a) G o~(a), the spectrum of (the 
operator) a for all a £ A. Moreover, |y?(a)| < ||ct|| and (p(a*a) > 0. 

• (Gelfand-Naimark theorem) Let A be an abelian C*-algebra and V the set of characters 
of A equipped with the weak* topology inherited from the dual A*, of A. It follows 
that V is a locally compact Hausdorff space which is compact if and only if A contains 
the identity. Moreover, A is isomorphic to the algebra C (V) of continuous functions 
over V which vanish at infinity. 

• A (symmetric) derivation 5 is a linear map from a *-subalgebra T>(5), the domain of 
S, into A with the properties that 

1. 5{A)* =5{A*), AeV(5) 

2. 8{AB) = 5{A)B + A5{B), A,B G V{5) 

• Theorem. ( \TdJ) Let A be a C*-algebra with identity 1, and let 5 be a norm-densely 
defined norm-closed operator on A with domain T>(5). It follows that S is the generator 
of strongly continuous one-parameter group of *-automorphism of A if and only if T>(5) 
is a *-algebra and 5 is a symmetric derivation. 



The listed properties of S imply 1 G T>(5) and 5(1) = 0. 
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• ( ||15|| ) Let again A = C(Q) be an abelian C*-algebra. By the dimension of Q we shall 
mean the Hausdorff dimension of Q. The classification results for derivations of A can 
be given in the following way: 

1. In the case that dimfl = it turns out that all closed derivations are trivial (i.e. 
equal to 0) if Q is total disconnected. 

2. The case dimfl — 1, A — C([0, 1]) has the complete classification. 

3. For dimension more than 2 only sporadic results are known. 

These results clearly show that in our discussion of the Hamilton picture (cf. Section 
II) implicitly we have used some additional topological properties of the phase space 

r. 

• Let Ai, i = 1, k and B be C*-algebras and Ai x • • -Ak the Cartesian product of 
Ai, Ak- A map $ : A\ x • • • x Ak — > B is said to be completely positive if, for every 
n > 0, the n-square £>-valued matrix [$(ai ;Pir , afc ;Pir .)]p r=1 is positive whenever the 
n-square ^-valued matrices [a»;p, r ]p ir are positive for % = 1, k. 

• Theorem ( j^J) Let $ : A — > B be a completely positive (in general, nonlinear) map 
between two C*-algebras A and B. Then, there exist (uniquely) completely positive 
maps $ m) „ : A — > B (m,n = 0,1,2,...) such that 

1. $(a) = Z)m=o E^=o ^m,n(fl) where a & A and the series in norm convergence. 

2- $ TO>n (za) = z m ^^ mn [a) where zG (C and a G i. The map $ m , n (-) with the 
above property is called (m, n)-mixed homogeneous. 

• Theorem ( j^j) Let $ : A — > B be (m, n)-mixed homogeneous map with m + n > 0. 
Then, there is a completely positive map <f> : A\ x • • • x ^4 m+ „ — >■ i3 with ^4 fe = ^. 

= 1, m + n such that 

$(a) = 0(a, a) (90) 

where aGi and 0(ai, a m + n ) is multilinear in (oi, a m ) and multi-conjugate-linear 
in (a m+ i, a m+n ). 
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